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Abstract 

The van der Waals and Casimir-Polder interaction of different atoms with graphene is investi- 
gated using the Dirac model which assumes that the energy of quasiparticles is linear with respect 
to the momentum. The obtained results for the van der Waals coefficients of hydrogen atoms and 
molecules and atoms of metastable He* and Na as a function of separation are compared with 
respective results found using the hydrodynamic model of graphene. It is shown that, regardless of 
the value of the gap parameter, the Dirac model leads to much smaller values of the van der Waals 
coefficients than the hydrodynamic model. The experiment on quantum reflection of metastable 
He* and Na atoms on graphene is proposed which is capable to discriminate between the two 
models of the electronic structure of graphene. In this respect the parameters of the phenomeno- 
logical potential for both these atoms interacting with graphene described by different models are 
determined. 

PACS numbers: 73.22.Pr, 78.67.-n, 12.20.Ds 
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I. INTRODUCTION 



It is common knowledge that dispersion force between two neutral atoms or molecules 
arises in second order perturbation theory from the dipole-dipole interaction This force 
takes its name from the fact that it is caused by the dispersions of dipole operators, i.e., 
by quantum fluctuations. At short separations, where the electromagnetic interaction can 
be considered as instantaneous, dispersion force depends only on one fundamental constant, 
the Planck constant H. In this case it is referred to as the van der Waals force Q]. At larger 
separations, where the retardation of the electromagnetic interaction becomes significant, 
dispersion force depends on both h and the velocity of light c. In this case the force is usually 
called the Casimir- Polder force 3| . As a result of interatomic interactions, dispersion forces 
act also between an atom (molecule) and a macroscopic body and between two closely spaced 
macroscopic bodies (in the retarded regime, the force between two macrobodies is called the 

n 

Casimir force |4(). 

The van der Waals and Casimir-Polder forces between atoms (molecules) and material 
waUs p.ay an important role in different physical, chemica. and biological phenomena fl 
|2j. The magnitude and the distance dependence of these forces were measured in several 
experiments j|J. In the last few years atom- wall interaction is being studied intensively in 
connection with the phenomenon of quantum-reflection |d-13|. Theoretical description of 

nmrn 

atom-wall interaction is given by the Lifshitz theory [l|, |2j, |4|, |5|, [14|] which expresses the van 
der Waals and Casimir-Polder energy and force through the dynamic electric polarizability 
of an atom and the frequency-dependent dielectric permittivity of wall material. A large 



body of research of atom- wall interactions for different atoms and wal 
of the Lifshitz theory has been performed in the literature [l, 2], 



ay 



materials on the basis 
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2l|. 



Considerable recent attention was attracted to two-dimensional carbon nanostructures, 
such as graphene, carbon nanotubes and fullerenes |22|. Carbon nanostructures possess 
unique mechanical, electrical and optical properties which are of major fundamental and 
applied interest. One of potential applications of carbon nanostructures is to the problem 
of hydrogen storage [23J . Keeping in mind that one-atom-thick nanostructures are not char- 
acterized by the dielectric permittivity, the Lifshitz theory seems to be not immediately 
applicable. Because of this, dispersion interaction between hydrogen atoms (molecules) and 
a graphite sheet (graphene) or single-walled carbon nanotubes was investigated mainly by 
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using the phenomenological density functional theory 24j-|27l|. The second-order perturba- 



tion theory was also used 



28j | to calculate line shifts of a two-level atom interacting with a 



nanotube. For the multiwalled carbon nanotubes with at least several walls, the concept of 



29]. 



dielectric permittivity of graphite was shown to be applicable 

The application of the Lifshitz theory of dispersion interaction to one-atom-thick carbon 
nanostructures calls for introduction of some other quantity which determines the reflection 
coefficients instead of usually used dielectric permittivity In order to extend the Lifshitz 
theory to the case of carbon systems, the description of graphene in terms of the two- 



dimensional free electron gas 



30] was used. Graphene was considered as an infinitesimally 



thin positively charged flat sheet, earring a homogeneous fluid with some mass and neg- 
ative charge densities (the hydrodynamic model). The interaction of the electromagnetic 



oscillations with such a sheet was described by the special reflection coefficients [31|, |32|. 
The proposed approximate description was applied for the calculation of dispersion interac- 
tion between two parallel graphene sheets [33], between graphene and material plate [341 ] . 
graphene and an atom or a molecule, and between an atom or a molecule and a single- walled 
carbon nanotube 35]. The hydrodynamic model, however, does not take into account that 



low-energy excitations in graphene are massless Dirac fermions except 
move with a Fermi velocity rather than with the speed of light 
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'or the fact that they 



361 ] . As a consequence, 



at low energies the dispersion relation for quasiparticles in graphene is approximately linear 
with respect to the momentum measured relatively to the corners of the graphene Brillouin 
zone. Under the assumption that these properties, which are referred to as the Dirac model, 
hold at any energy, the reflection coefficients for the electromagnetic oscillations on graphene 
were found 0- These coefficient s were used to calculate the Casimir interaction between a 
graphene and a parallel ideal metal plane {37 ]. 

In this paper we investigate the dispersion interaction of atoms (molecules) with graphene 
described by the Dirac model. We consider an atom and a molecule of hydrogen, and also 
atoms of metastable He* and Na often used in experiments on quantum reflection. The van 
der Waals coefficient C3 is computed by using the Lifshitz formula as a function of atom- 
graphene separation. The obtained results are compared with those obtained previously 
using the hydrodynamic model of the electronic structure of graphene. It is shown that in 
all cases the magnitudes of the van der Waals and Casimir- Polder energy computed with the 
help of the Dirac model are less than with the hydrodynamic model. Differences between the 
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two sets of results increase with the increase of separation. We also find the dependence of 
the predicted van der Waals and Casimir- Polder energies on the value of the gap parameter. 
We propose the experiment allowing to discriminate between the predictions of the Dirac and 
hydrodynamic models of graphene. For this purpose the parameters of phenomenological 
potential of atom-graphene interaction are determined for atoms of metastable He* and Na. 

The paper is organized as follows. In Sec. II we briefly present the Lifshitz formula for 
atom-wall interaction and the forms of reflection coefficients in both models. Section III 
contains computational results for the van der Waals coefficient C3 as a function of sepa- 
ration and gap parameter for atoms and molecules of hydrogen. In Sec. IV similar results 
are presented for atoms of metastable He* and Na. The parameters of phenomenological 
potential for these atoms are also determined for the needs of proposed experiment. In 
Sec. V the reader will find our conclusions and discussion. 



II. TWO MODELS FOR INTERACTION OF ATOMS WITH GRAPHENE 



The van der Waals and Casimir- Polder energy in the configuration of a microparticle (an 
atom or a molecule) and a plane structure can be expressed in terms of reflection coefficients 
of the electromagnetic oscillations on this structure in the following way 4], 



14 



E(a) = / <i£a(i£) / k±dk±q< 

2vr Jo * Jo 



-2ag 



x <^ 2r TM (i£, k x ) ^ [ruiK, k ±) + r TE(^, k±)] 



q 2 c 2 



(1) 



Here, a is the separation distance between the microparticle and the plane, k± is the mag- 
nitude of the wave vector projection on the plane (the latter is perpendicular to the z 
axis), a(i£) is the dynamic polarizability of a microparticle calculated along the imaginary 
frequency axis, and q 2 = k ± + $, 2 /c 2 . The reflection coefficients for two independent po- 
larizations of the electromagnetic field, transverse magnetic and transverse electric, ?"tm,tE) 



are also calculated along the ima 
approach (see, for instance, Refs 



bi nary frequency axis. In the framework of the scattering 
4J,l38[]) it can be proved that under some conditions Eq. (CQ) 



holds for any planar structure with appropriately chosen reflection coefficients. Specifically, 
the use of Eq. (CQ) is justified if a plane structure (a material of the plate or a graphene) is in 
thermal equilibrium with an environment at some not too high temperature, and an atom is 
in the ground state [l5(. For instance, the probabilities that an atom is in an excited state 



are distributed according to the Boltzmann law and are negligibly small at room temper- 
ature (T = 300 K). As a result, the zero-temperature Lifshitz formula ([1]) can be used at 
sufficiently small separations (see below for quantitative evaluation of its application region). 

There are different models of reflection coefficients for graphene using different boundary 
conditions. In the framework of the hydro dynamic model of the electronic structure of 
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35] 



graphene they can be presented in the form 



c 2 qK 



± 



c 2 qK + i 2 ' 
K 



(2) 



K + q 

where the wave number of the graphene sheet K = 6.75 x 10 5 m _1 corresponds to the 
frequency uk = cK = 2.02 x 10 14 rad/s. As was explained in Sec. I, the hydrodynamic 
model does not take into account some important properties of graphene which hold at low 
energies, specifically that the energy of quasiparticles is a linear function of the momentum. 
These properties are well described by the Dirac model of the electronic structure of graphene 
where the reflection coefficients are given by 37] 



aq&(q) 



2q 2 + aq&(q)' 
a®(q) 



(3) 



2q + a$(q) ' 

Here a = e 2 /(he) ~ 1/137 is the fine-structure constant, q 2 = (v^k^ + £ 2 )/c 2 , v-p ~ 10 6 m/s 
is the Fermi velocity, and the function $ determines the polarization tensor in an external 
electromagnetic field in the one-loop approximation in three dimensional space-time. The 
explicit form of this function along the imaginary frequency axis is the following 37] 



= N A + 



-4A 2 q 

arctan 

2q 2A 



(4) 



where for graphene N = 4, A = A/ (he), and the value of the gap parameter A is not well 
known. The upper bound on A is approximately equal to 0.1 eV, but it might be also much 



smaller 
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Below, we perform computations of the van der Waals and Casimir-Polder interaction 
energy between different atoms and graphene. For this purpose it is useful to introduce the 
dimensionless variable y = 2qa instead of k±. Then Eq. ([1]) can be identically rewritten in 
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the form 



E(a) 



CM 



(5) 



where the van der Waals coefficient Cs(a) is given by 



C*(a) 



h 

167T 



dye' 



cy/(2a) 



x <^ 2y 2 r TM (i£,y) 



2c2 



[rn*(i£,y)+rTE(i€,y)] 



(6) 



Either the hydro dynamic- or Dirac-model reflection coefficients (j2]) and ([3]), respectively, can 
be substituted here keeping in mind that in terms of new variable it holds 



Q 



y_ 

2a 



q 



y 2 

c 2 Aa? 



+ 1- 



e 



1/2 



(7) 



Notice that Eq. fl5]) with constant C3 presents the energy of nonretarded van der Waals 
interaction between an atom and a wall which is valid only at small separations of about 
nasi , ll4H16| (at shorter separations the description of graphene by means of 



2-3 nm [ 

the boundary conditions becomes inapplicable and one should take into account an atomic 
structure of the sheet). If, however, C3 depends on separation in accordance with Eq. (jBJ), the 
interaction energy (jSJ) includes both the nonretarded and retarded (relativistic) regimes. As 
a result, Eqs. (JSJ) and (JSJ) are applicable up to separation distances where the thermal effects 
become significant. It is well known J4, 15 1 that at room temperature the zero-temperature 
Lifshitz formula provides very exact computational results of atom-plate interaction up to 
separations of about 1 /xm (see Sec. IV for a computational example). As to the function 
$ determining the polarization tensor, relativistic thermal quantum field theory does not 
predict a noticeable change of this quantity at room temperature, as compared to the case 
of zero temperature (note that computation in Ref. [39| was, strictly speaking, performed 
in a nonretarded regime). 



III. INTERACTION OF HYDROGEN ATOMS AND MOLECULES WITH 
GRAPHENE 

Here, we compare the van der Waals coefficients for interaction of hydrogen atoms and 
molecules with graphene in the cases when computations are performed with the help of 
Dirac and hydrodynamic models. To perform computations using Eq. (jHJ), one needs some 

6 



expressions for the atomic and molecular dynamic polarizabilities of hydrogen. As was 



shown in Refs. 
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29], for the calculation of dispersion forces, the polarizabilities can be 



represented with sufficient precision in the framework of single-oscillator model 

.... a(0) 

«(*0 = — 4r, ( 8 ) 

where a(0) is the static polarizability and uo is the characteristic frequency. For a hydrogen 
atom it was found [40] that a(0) = a a (0) = 4.50 a. u. and u>o = u>o a = 11.65 eV (remind that 

. a.u. of polarizability is equal to 1.482 x 10~ 31 m 3 ). For a hydrogen molecule it holds that 

4(J a(0) = a m (0) = 5.439 a.u. and oo = u 0m = 14.09 eV. 

In Fig. 1(a) we present the computational results for the van der Waals coefficient C33 
(in atomic units) for a hydrogen atom interacting with graphene as a function of separation. 
Note that one atomic unit for C3 is equal to 4.032 x 10~ 3 eVnm 3 . Computations were 
performed using Eqs. ([6]) and (jSJ) with the hydrodynamic-model reflection coefficients (jSJ) (the 
dashed line) and with the Dirac-model reflection coefficients ([3]) (the upper and lower solid 
lines obtained with the lower and upper bounds for graphene gap parameter A, respectively). 
As an upper bound, the value 0.1 eV was chosen (see Sec. II). Keeping in mind that below 
some A min the values of C3 are insensitive to further decrese of A (see below), the value of 
10 _15 eV was chosen as a lower bound. As is seen in Fig. 1(a), the values of C 3j h computed 
using the hydrodynamic and Dirac models are significantly different and the ratio C^/C^ 
increases with the increase of separation distance. Thus, at the shortest separation a = 3 nm 
it holds C^fi/C^fi = 1.065. When separation varies to 5, 10, 20, 50, and 100 nm this ratio 
increases to 1.19, 1.44, 1.85, 2.85, and 4.21, respectively (we have used the computational 
data related to the lower solid line with A = 0.1 eV). At large separations of about 100 nm 
both Cgg(a) and C^(a) decrease with separation as 1/a leading to E# (a) ~ 1/a 4 . This 



is a typical behavior of the atom-plate interaction at relativistic separations 3)- 5] which 
holds also for other atomic systems considered below. In the nonretarded region a < 3nm 
the Casimir-Polder energy depends on separation as or 7 / 2 . This region is, however, outside 
the application region of our formalism based on the use of boundary conditions. 

The relationship between the two solid lines in Fig. 1(a) demonstrates that the influence 
of the gap parameter on the value of C 3j h is rather moderate. This is illustrated in more 
detail in Fig. 1(b) where C 3j h is plotted as a function of log 10 A. The solid lines labeled 
1, 2, and 3 are computed at separations a = 5, 50, and 100 nm, respectively. For the 



line labeled 1 the relative difference between the maximum and minimum values of C 3i h, 
(C££ x - C™£)/C™£ = 6.6%. For the lines labeled 2 and 3 this relative difference is equal to 
16.4% and 31.3%, respectively. For sufficiently small A min (different at each separation) C 3j h 
becomes constant and does not depend on further decrease of A. Thus, for lines labeled 1, 2, 
and 3 the van der Waals coefficient is constant for gap parameters satisfying the inequalities 
A < 0.01 eV, A < 0.004 eV, and A < 0.001 eV, respectively. 

The computational results for a hydrogen molecule differ from the respective results 
for a hydrogen atom quantitatively but not qualitatively. In Fig. 2 the van der Waals 
coefficient C3 h 3 is plotted as a function of separation [the dashed line is computed using 
the hydrodynamic model and the two solid lines using the Dirac model with the same gap 
parameters as in Fig. 1(a)]. Here, at the shortest separation a = 3nm it holds C3H /^3H 2 = 
1.045, i.e., the computational results are slightly closer than for a hydrogen atom. With the 
increase of separation to 5, 10, 20, 50, and 100 nm the ratio of the van der Waals coefficients 
computed using the two models increases to 1.18, 1.45, 1.89, 3.00, and 4.63, respectively, 
i.e., becomes larger than respective ratio for a hydrogen atom. As to the role of the gap 
parameter of graphene, its value only slightly influences the computational results. Figures 
1(a) and 2 raise a question if it is possible to discriminate between the two models of graphene 
from experiments on quantum reflection. This question is considered in the next section. 



IV. QUANTUM REFLECTION OF METASTABLE He* AND Na ATOMS ON 
GRAPHENE 



Experiments on quantum reflection 



4l3| of cold atoms on different surfaces demon- 



strated that the reflection amplitude depends critically on the form of interaction potential. 
Because of this, using the experimental data for reflection amplitudes, it is possible to re- 
store the form of potential within the limits of experimental error. If the two models of 
graphene (the hydrodynamic one and the Dirac one) lead to the values of interaction energy 
which differ dramatically, it might be possible to experimentally exclude some model or find 
it consistent with the data. To make reflection amplitudes more pronounced it is customary 
in experiments on quantum reflection to use atoms with large static electric polarizability. 
This allows to obtain more than one order of magnitude larger values of the van der Waals 
coefficient C3 when compared to the values of the same coefficient for hydrogen atoms and 
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molecules plotted in Figs. l(a,b) and 2. 

We begin with an atom of metastable He* interacting with graphene by means of 
dispersion interaction. For He* atom it holds [4l| a(0) = an c *(0) = 315.63 a. u. and 
Uq = w ,He* = 1.18 eV. Computations were performed using Eqs. (EJ) and ([HD with reflec- 
tion coefficients, as defined in the hydrodynamic 02]) or Dirac fl3]) models of graphene. The 
computational results for the van der Waals coefficients C3. 

He* cLS cL function of separation are 
presented in Fig. 3 by the dashed line (the hydrodynamic model) and by the two solid lines 
(the Dirac model) obtained with two values of the gap parameter of graphene indicated in 
Sec. III. Similar to the case of hydrogen atom and molecule the ratio C^ e */C^ e * ls larger 
than unity (it is equal to 1.33 at a = 3nm and increases with the increase of separation). 
At a = 5, 10, 20, 50, and 100 nm this ratio is equal to 1.47, 1.76, 2.23, 3.33, and 4.78. 
Thus, it is larger than for atoms and molecules of hydrogen at all separations. What is 
more important, the magnitudes of C^ e * computed with the Dirac model are by a factor 
of 18 (at a = 3 nm) and by a factor of 46 (at a = 100 nm) larger than respective values of 
Cgg computed using the same model. In a similar way, C^^^nm) = 23(7.^(3 nm) and 
Cg^ e *(100nm) = 52C3^(100 nm). This makes experiments on quantum reflection of He* 
atoms on graphite feasible. Note that the role of thermal effects at separations considered 
is negligibly small. For example, we have computed the values of C^ c „ at a = 0.5 fim 
both at T = and T = 300 K using the zero-temperature and the thermal Lifshitz for- 
mulas. In the framework of the Dirac model the following respective values were obtained: 
^3 He* = 0.0183505 a.u. and C^ e * = 0.0183565 a.u. This means that at a = 0.5 fim the rela- 
tive difference between the results obtained using the zero-temperature and thermal Lifshitz 
formula is equal to only 0.033%. 

It is interesting also to analyze the influence of a gap parameter of graphene on the value 
of C 3i He* illustrated in Fig. 3 by the upper and lower solid lines computed with A = 10~ 15 eV 
and with A = 0.1 eV, respectively. Thus, at a = 5nm it holds (C 3 m ^, - C™™ e *)/C™™ e * = 
3.7%. At separations a = 50 and 100 nm the same relative difference is equal to 28% and 
46%, respectively. From this it follows that at large separations the impact of the value of 
a gap parameter A on C^ e *(a) is somewhat stronger than on C^(a). 

Another atom suitable in experiments on quantum reflection on graphene is Na. For Na 



atom it holds |42| a(0) = aN a (0) = 162.68 a.u. and u = u; 0i Na = 1-55 eV. The computational 
results for the van der Waals coefficients C^ a (the dashed line) and C^ a (the two solid 
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lines for two values of the gap parameter indicated above) as a function of separation are 
presented in Fig. 4. At a = 3nm we obtain Cf^/C^l = L4 °- With the 

increase of 

separation to a = 5, 10, 20, 50, and 100 nm this ratio increases to 1.55, 1.87, 2.39, 3.61, 
and 5.29, respectively. The magnitudes of are by a factor 11.6 (at a = 3nm) and by 
a factor of 26 (at a = 100 nm) larger than the respective values of Cg^. In a similar way, 
C^ a (3nm) = 150^(3 nm) and cg} a (100nm) = 330^(100 nm). The exceeding of C^ a 
relative to Cgg is not so large as for He* atoms, but still sufficient for experimental purposes. 

The upper and lower solid lines in Fig. 4 illustrate the influence of a gap parameter of 
graphene A on the value of C 3iNa . Thus, at a = 5nm it holds (C 3 m ^ - Cg^/Cgf^ = 3.2%. 
At separations a = 50 and 100 nm the relative difference between maximum and minimum 
values of C^Na is equal to 25.6% and 42%, respectively. These results are similar to the 
respective results for He* atoms. 

In almost all papers devoted to the experimental investigation of quantum reflection 



0,E 



ii. 



13] the comparison of the measurement data with theory is performed with the use 



of a phenomenological potential for the atom-wall interaction. This simplifies the calculation 
of theoretical reflection amplitudes to be compared with the exp erimental ones. The most 
often used phenomenological potential has the form (4 -7, 
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13] 



where I is a characteristic parameter with the dimension of length that depends on the 
material of the plate. It is assumed that at short separations a <C / (typically at separations 
of a few nanometers where the retardation effects are negligibly small), E ph (a) coincides with 
the van der Waals interaction energy between an atom and usual material plate [Eq. (|5]) with 
practically constant C3], so that C4 = IC3. At large separations, where I is negligibly small 
in comparison with a, the phenomenological potential OH]) coincides with the Casimir-Polder 
retarded interaction energy — C^/a 4 obtained for ideal metal wall The comparison 

of computational results obtained using Eq. (Q with measurement data for the reflection 
amplitudes allows one to determine the experimental values for the parameters / and C4. 
Note that Eq. does not agree with a nonretarded asymptotic limit for atom-graphene 
interaction (see Sec. III). As shown below, Eq. (|9]) provides a very good agreement with 
computations using formula o, Jy at" oderate Rations. In Re, Q ^ 

accuracy of the phenomenological potential was investigated by the comparison with exact 
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results obtained using the Lifshitz theory. It was shown that the phenomenological potential 
leads to quantitatively correct results only in the limit of short separation distances below 
100 nm. This conclusion was confirmed experimentally in Ref. 43] where deviations between 
the theoretical prediction using the phenomenological potential and the measurement data 
within the separation region from 160 to 220 nm were demonstrated. Below we use the 
phenomenological potential (Q only at separations below 100 nm and compare the obtained 
results with those using the Lifshitz theory. The separation regions where Eq. (Q leads to 
less than 1% error are specified. 

Here, we perform the best fit of the phenomenological potential (Q to the energies E^ h \a) 
and E^ D \a) computed above for the interaction of He* and Na atoms with graphene de- 
scribed using the hydrodynamic and Dirac models, respectively. In doing so we first use the 
gap parameter of graphene A = 0.1 eV. The results of the best fit of Eq. (Q to the dashed 
line in Fig. 3 (a He* atom interacting with graphene described by the hydrodynamic model) 
are the following: 

C i% e * = 85.11 a.u., zjjj. = 72.77 nm. (10) 

We note that 1 a.u. of C4 is equal to 4.032 x 10 -3 eVnm 4 . For the Dirac model, by making 
the best fit of Eq. (Q to the lower solid line in Fig. 3, we obtain 

Cg£. = 12.59 a.u, = 11.18 nm. (11) 

It can be seen that the value of E ph (a) from Eq. ([9]) with the parameters ( flOj) deviate from 
the dashed line in Fig. 3 for less than 1% in the separation region from 10 to 60 nm. At the 
shortest and largest separations (3 and 100 nm, respectively) relative deviations achieve 10%. 
For the Dirac model, Eq. (Q with the parameters ffTTl) deviates from the lower solid line in 
Fig. 3 for less than 1% at separations from 6 to 100 nm. The largest relative deviation of 
5% holds at a = 3 nm. So large differences in the parameters of phenomenological potential, 
as predicted by the two models of graphene (by a factor of 6.8 for CiHe* an d by a factor 
of 6.5 for Zne*)) make possible to decide between these models experimentally. The use of 
alternative values of the gap parameter leaves this conclusion unchanged. For example, 
repeating the same fit with A = 10 _15 eV (the upper solid line in Fig. 3) one arrives at 
^4 He* = 18.04 a.u. and l^J, = 18.22 nm. These are by factors 4.7 and 4.0 smaller than 
respective values (fTOl) obtained using the hydrodynamic model of graphene. 
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Similar fit of Eq. (Q was performed to the lines of Fig. 4 representing the computational 
results for atoms of Na. This leads to 

Cg> a = 50.82 a.u., ZjJi = 66 - 92 nm ( 12 ) 

for the hydrodynamic model and 

C$ a = 7.11 a.u., 4? = 9.77nm (13) 

for the Dirac model of graphene. The energy (Q with parameters ([12]) deviates from the 
dashed line in Fig. 4 for less than 1% in the separation region from 10 to 60 nm. At 
separations of 3 and 100 nm the relative deviations achieve 11%. For the Dirac model, 
Eqs. ([H]) and (1131) describe the lower solid line in Fig. 4 with less than 1% error in the region 
from 6 to 100 nm with maximum relative deviation of 5.7% at a = 3nm. Thus, for Na 
atoms the predictions of the two models of graphene for C 4j Na and differ by factors of 
7.1 and 6.8, respectively. This makes possible an experimental choice between these two 
models. The above computations were performed with A = 0.1 eV. For the extremely small 
gap parameter A = 10~ 15 eV (the upper solid line in Fig. 4), one arrives at the values 
^4Na = 9-74 a.u. and l^J = 15.45 nm. These parameters of the phenomenological potential 
09]) are still less than the parameters obtained with the use of hydrodynamic model by factors 
of 5.2 and 4.3, respectively. This confirms that the experimental choice between the two 
models is possible with any value of the gap parameter. 



V. CONCLUSIONS AND DISCUSSION 



In the foregoing, we have investigated the van der Waals and Casimir-Polder interaction 
of different atoms with graphene on the basis of fundamental Lifshitz theory of atom-wall 
interaction. This was done by using the two approximate models of the electronic structure 
of graphene proposed in the literature, the hydrodynamic model and the Dirac model. The 
respective two sets of reflection coefficients for the electromagnetic oscillations on graphene 
were used in computations of the van der Waals coefficients describing dispersion interaction 
of hydrogen (H and H2), and He* and Na atoms with graphene as a function of separation. 
It was shown that the hydrodynamic and Dirac models of graphene lead to drastically 
different predictions for the strength of atom-graphene interaction. This conclusion was 
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obtained independently of the value of the gap parameter of graphene which was allowed to 
reduce from the experimental upper limit to almost zero. 

Keeping in mind that both models of graphene are idealizations valid for different values 
of parameters, the new experiment was proposed capable to discriminate between the alter- 
native theoretical predictions. For this purpose it was suggested to consider the quantum 
reflection of He* and Na atoms on graphene. Due to large static electric polarizabilities, 
the van der Waals coefficients for the interaction of He* and Na atoms with graphene were 
shown to be in several tens times greater than for H and H 2 . This makes the observation 
of quantum reflection of He* and Na atoms on graphene feasible. For the needs of proposed 
experiment, we have determined the parameters of phenomenological potential for He* and 
Na atoms interacting with graphene described by both the hydrodynamic and Dirac models. 
It was shown that the parameters obtained using the hydrodynamic model are by a factor 
varying from 4.0 to 7.1 greater than the parameters obtained using the Dirac model. This 
demonstrates that the proposed experiment should be capable to discriminate between the 
predictions of the Dirac and hydrodynamic models of graphene. 

In the future it is supposed to employ the Dirac model for the computation of dispersion 
interaction between graphene and material plates made of different materials and between 
various atoms and molecules and carbon nanotubes. 
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FIG. 1: The van der Waals coefficient for a H atom interacting with graphene as a function of (a) 
separation and (b) gap parameter. The dashed and the two solid lines are computed using the 
hydrodynamic model and the Dirac model with two different values of the gap parameter indicated 
in the text. The solid lines labeled 1, 2, and 3 are computed using the Dirac model at separations 
a = 5, 50, and 100 nm, respectively. 
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FIG. 2: The van der Waals coefficient for a H2 molecule interacting with graphene as a function of 
separation. The dashed and the two solid lines are computed using the hydrodynamic model and 
the Dirac model with two different values of the gap parameter indicated in the text. 
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FIG. 3: The van der Waals coefficient for a He* atom interacting with graphene as a function of 
separation. The dashed and the two solid lines are computed using the hydrodynamic model and 
the Dirac model with two different values of the gap parameter indicated in the text. 
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FIG. 4: The van der Waals coefficient for a Na atom interacting with graphene as a function of 
separation. The dashed and the two solid lines are computed using the hydro dynamic model and 
the Dirac model with two different values of the gap parameter indicated in the text. 
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